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(54) Methode pour generer un maillage 3D 

(57) -La methode peut s'appliquer par exemple a 
uno structure geologique definie a partir de don- 
noos disponibles acquises soit par des mesures, 
soit par des analyses et des interpretations de don- 
nees sismiques ou des diagraphies realisees dans 
des puits. Elie comporte une modelisation prealable 
des limites (bords de la zone, discontinuity ou 
failles par exemple) par des surfaces triangulees, 
enveloppant au mieux ces limites, et un fractionne- 
ment de la zone d'etude en blocs ayant au plus six 
faces, dont les faces correspondent a des disconti- 
nuites ou des failles principales de la zone preala- 
blement triangulees. Chaque bloc est ensuite sub- 
divise en un ensemble de cellules elementaires ne- 
cessaires pour la modelisation de la zone par ele- 
ments ou differences finis, en realisant un maillage 



de chacune des faces du bloc par interpolation a 
partir de ses bords defacon a obtenir une grille, sui- 
vi d'un peuplement en noeuds de chaque bloc dans 
son epaisseur, par interpolation a partir des noeuds 
des diff erentes grilles formees sur les diflerentes fa- 
ces, ces operations etant accompagnees de proce- 
dures particulieres d'homogeneisation pour reajus- 
ter la taille des cellules a I'interieur de chaque bloc 
suivant des regies definies. Les blocs mailles sont 
ensuite reassembles. 

Application notamment a I'etude de gisements pe- 
troliferes. 
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Description 

La presente invention concerne une methode pour generer un maillage 3D dans le but de modeliser un corps 
quelconque dont la geometrie est connue et imposee a priori tel que par exemple une structure geologique en tenant 
5 compte fidelement de sa geometrie dans toute sa complexity naturelle. Le maillage realise par la methode s'accorde 
bien dans le cas d'une structure geologique, avec la position des interlaces entre les couches, avec les failles meme 
de forme complexe, les biseaux et autres irregularites de structure, en les representant avec leur conformation reelle. 

La methode selon I'invention trouve des applications notamment pour la modelisation en geomecanique : pour la 
caracterisation des reservoirs : couches de sediments, stratigraphie, et aussi pour la simulation des ecoulements dans 
to les bassins et reservoirs. 

I) ETAT DE LA TECHNIQUE 

La modelisation d'une structure geologique est realisee a partir d'un ensemble de donnees disponibles acquises 

is soit par des mesures, soit par des analyses et des interpretations de donnees sismiques ou des diagraphies realisees 
dans des puits. Ces donnees se presentent sous la forme d'un nuage de points entre lesquels, au mieux, il convient 
d'installer un maillage. La formation d'un maillage en 3D joue un role crucial pour faire des simulations numeriques 
dans differentes applications geologiques. Beaucoup de recherches ont ete menees sur le sujet principalement en 
ingenierie de reservoir en raison de ses applications dans le domaine petrolier. 

20 On connait de nombreuses techniques de maillage 3D pour des simulations numeriques par elements ou diffe- 

rences finis et des ensembles de logiciels de simulation crees pour leur mise en oeuvre qui sont adaptes a des appli- , 
cations et respectent des contraintes determinees. 

II est connu pour certaines applications d'utiliser des maillages de type tetraedrique non structure. II est connu 
aussi de faire du maillage hexaedrique de. peau ou de surface et de repartir des points dans le maillage sort par 

25 interpolation transfinie, soit encore par resolution d'equations aux derivees partielles. Les surfaces sont generalement 
creees a partir de fonctions parametriques et les points du maillage sont obtenus simplement par echantillonnage de 
ces fonctions. Une telle methode ne convient pas generalement pour des applications dans les geosciences ou les 
volumes a modeliseront des formes le plus souvent tres complexes, difficiles a traduire par des fonctions parametriques 
compte-tenu de la precision qui est recherchee Aussi prefere ton approcher les formes complexes par des elements 

30 de surface triangules. 

Differentes techniques de maillage appliquees a la modelisation de reservoir, sont decrites par exemple dans : 

Georges P.L : "Generation Automatique de Maillages", Editions Masson, 1991; 

35 - Guerillot D R. et al: An ineractive 3D Mesh Builder for Fluid Flow Reservoir Simulation: in SPE - Computer Appli- 
cations, Dec. 1993; 

Aziz K.; Reservoir Simulation Grids : Opportunities and problems. SPE-JPT, July 1993; ou 

40 - Heineman Z.E. et al : Gridding Techniques in Reservoir Simulation. Intl Forum on Reservoir Simulation. Alpbach 
Austria, Sept. 12-16, 1983. 

Cependant, la plupart des techniques de maillage connues, et de la les simulations qui en resultent, ne conviennent 
bien que pour modeliser des structures geologiques relativement simples. Les irregularites geologiques courantes 

45 dans les zones-reservoir modelisees, telles que les couches ayant de grandes variations d'epaisseur, les failles ou les 
biseaux.. sont rarement pris en compte ou alors seulement sous une forme grandement simplifiee. Dans les modeles 
mailles obtenus dans le cadre des geosciences, les failles apparaissent verticalement; au mieux, elles sont inclinees 
el planes. Or il est important de rendre fidelement, la geometrie en 3D des sediments car ils possedent des caracte- 
ristiques physiques propres, et aussi la geometrie des failles et autres discontinues mecaniques. 

so La methode selon I'invention permet de generer un maillage 3D respectant la geometrie d'une zone d'un milieu 

complexe dont la geometrie est connue et imposee a priori, tel qu'une structure geologique, dans le but de realtser un 
modele representatif de ce milieu. La structure est .definie a partir de donnees disponibles acquises soit par des me- 
sures, soit par des analyses et des interpretations de donnees sismiques ou des diagraphies realisees dans des puits. 
Elle comporte une modelisation prealable des limites (bordsde la zone, discontinuites ou failles par exemple) par des 

55 surfaces triangulees, enveloppant au mieux ces limites, et un fractionnement de la zone d'etude suivant en blocs ayant 
au plus six faces, dont les faces correspondent a des discontinuites ou des failles principales de la zone prealablement 
triangulees. 

Elle est caracterisee en ce que'elle comporte : 
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la subdivision de chaque bloc en un ensemble de cellules elementaires necessaires pour mpdeliser la zone par 
elements finis, par une reticulation en 3D de chacun des differents blocs obtenus, subdivision qui comprend : 

a) le maillage de chacune des faces du bloc par interpolation a partir de ses bords, de facon a creer une grille 
intermediaire, suivi d'une projection sur la surface triangulee des noeuds de cette grille intermediaire suivant 
une direction definie (perpendiculairement a elle-meme par exemple); 

b) le peuplement en noeuds de grille de chaque bloc dans son epaisseur, par une methode d'interpolation a 
partir des noeuds des differentes grilles formees sur les differentes faces triangulees de chaque bloc, et 

c) un rearrangement dans le positionnement des noeuds des cellules elementaires creees, obtenu (de facon 
iterative, generalement) par relaxation dans le positionnement des noeuds, en fonction de criteres predefinis 
tels que la regularity le choix de la taitle relative des differentes cellules et le mode de progression de leurs 
tallies, la direction de projection etc. 



Chaque bloc etant ainsi reticule en 3D, on peut proceder a un reassemblage ou reempilement de tous les blocs. 
On peut proceder successivement a une premiere operation de relaxation dans le positionnement des noeuds sur 
chacune des differentes faces triangulees, le peuplement en noeuds internes etant effectue a partir des grilles relaxees, 
et ensuite a une deuxieme operation de relaxation dans le positionnement de ces noeuds internes. 
20 De preference, la premiere operation de relaxation est effecluee par des phases ileratives de relaxation-projection. 

La methode selon I'invention permet d'obtenh un modele qui respecte ties precisement la geometrie souvent 
complexe de structures geologiques telles que des champs petroliferes, ce qui a pour effet que les simulations qui 
peuvent etre faites a partir de ce modeleVendent plus fidelement compte des phenomenes. La methode utilisee permet 
en outre une economie tres importante de temps de calcul pour realiser les operations de maillage 3D habituellement 
25 tres longucs avec les rnethodes existantes. 

D'autres caracteristiques et avantages de la methode selon I'invention, apparaTtront a la lecture de la description 
ci-apres de modes de realisation decrits a titre d'exemples non limitatifs, en se referant aux dessins annexes ou : 

' - les Fig. 1 et 2 montrent des exemples de relaxation-projection sur un maillage surfacique 4X4; 

30 

la Fig. 3 montre un cube deforme avec des faces triangulees; 

la Fig. 4 montre t'enveloppe externe d'un cube maille; 

35 - les Fig. 5 et 6 montrent les surfaces internes d'un cube maille respectivement avant et apres relaxation; 

la Fig. 7 montre des courbes relevees sur cinq coupes d'une structure geologique de tres grandes dimensions, sur 
laquelle la methode selon I'invention a ete testee, qui constituent une partie des donnees initiales dont on dispose, 
sur cette structure; 

40 

la Fig. 8 montre d'autres donnees initiales sur la meme structure, en I'occurrence, des cartes d'isobathes et les 
emplacements de 4 puits; 

la Fig.9 montre une modelisation surfacique ( a base de surfaces triangulees) de la structure geologique etudiee; 

45 

la Fig. 10 montre un exemple de decoupage de la structure en blocs; 

la Fig. 11 montre le resullal d'un maillaye bloc par bloc de I'enveioppe externe de la structure geologique; 

50 - [a Fig. 12 montre le resultat d'un maillage bloc par bloc des surfaces internes de la structure geologique; 

la Fig. 13 illustre a titre d'exemple comment on prend en compte un biseau dans le maillage de la structure geo- 
logique; 

55 - la Fig. 14 illustre a titre d'exemple comment on prend en compte une faille dans le maillage de la structure geolo- 
gique; 

la Fig. 15 illustre la methode d'interpolation 2D de Coons; 
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la Fig. 16 montre les elements pris en compte dans la methode d'interpolation 3D de Coons; et 

la Fig. 17 illustre le mode de relaxation-projection des noeuds de maille. 

5 La structure geologique de grandes dimensions illustree par les Fig. 7-10, correspond a un dome anticlinal com- 

portant six interfaces entre des sediments. 

La realisation d'un maillage en 3D par elements finis d'une zone a partir d'un ensemble de donnees initiales con- 
nues pour la zone : interfaces entre differentes couches ou horizons, failles principales, carte en courbes de niveau 
avec les emplacements de puits passe par une succession d'etapes qui seront decrttes ci-apres. 

io 

I) Moderation des limites ou bords de la zone etudiee et de ses principales discontinuites. 

On realise une modelisation des surfaces : bords, interlaces, failles en utilisant une logiciel specialise connu tel 
que GOCAD (marque deposee) sous la forme d'une juxtaposition continue d'elements de surface triangulares. Dans 
is un premier temps on forme une ebauche qui s'adapte au mieux aux donnees initiales. De maniere iterative, on deplace 
les sommets de la surface en utilisant un logiciel d'interpolation et on densitie localement les elements de surface 
triangulares, de facon a plaquer au mieux la surface sur le semis de points (donnees d'entree). Les elements de 
surface ainsi crees ne se raccordent generalement pas comme ils le devraient : les interfaces ne sont pas limitees aux 
intersections et au bord du domaine d'etude et elles peuvent se croiser en debordant. On calcule alors les lignss 
d^nterception entre surfaces et on limite celles-ci aux courbes d'interseclion. La Fig 7 monlre le resultat final de celle 
in_>delisation des limites de la structure et Ton peut remarquer que la plupart des triangles des surfaces triangulees du 
biicr:u sont localises dans les regions de forte courbure. 

A ce stade, il n'est pas encore possible de produire un maillage 3D tenant compte des contraintes geologiques et 
numcnques. Pour y parvenir, on procede a un fractionnement de la zone. 

II ) Fr -ictionnement de la zone etudiee 

On decoupe d'abord la zone en macro-blocs ayant au plus six faces. II s'agit generalement de blocs de fcrme 
hexaedrique (a six faces), mais le nombre de faces peut etre inferieur a six et descendre jusqu'a deux. Dans ce cas, 

30 pour les besoins de la methode, le bloc est quand meme considerere comme un block hexaedrique mais ayant des 
faces degenerees (reduites a un point ou a une courbe). Le mode de decoupage est generalement suggere par la 
geometrie des horizons, des principales failles et leur interaction. Pour constituer ces macro-blocs 3D, on a besoin 
parfois d'introduire des frontieresartificielles. C'est necessatre par exemple quand un bloc a plus de six faces naturelles 
et ne satisfait done pas a la regie choisie d'etre au plus hexahedrique. Dans un tel cas : on reunit deux faces en une 

55 face complexe pour satisfaire le critere, ou bien on introduit des surfaces de coupure artiftcielles. 

Un exemple de decoupage de la structure geologique etudiee en macro-blocs est montre a la Fig. 10, ou Ton a 
extrapole la gaille jusqu'a la frontiere de la structure et verticalement projete I'intersection du biseau sur les autres 
surfaces. 

40 ill) Maillage 

On realise une reticulation ou maillage en 3D separement de chacun de ces macro-blocs de fagon a le subdiviser 
en blocs elementaires qui vont former les cellules pour la modelisation par elements finis recherchees. Cette operation 
est realisee en deux temps. 

45 

1) Maillage des faces 

a) on cree tout d'abord une grille sur chacune des cinq ou six faces de chaque macro-bloc. Pour cela, on forme 
une grille intermediaire par interpolation a partir des quatre bords de la face, en utilisant la methode d'interpolation 
so transfinie de Coons dont le principe est rappele ci-apres. 

Etant donne -iuatre courbes parametriques fl(u), f2(u), gl(v), g2(v) (0 2 u, v 2 1) definissant quatre bords jointifs 
(Fig.15) : on introduit une forrnule pour calculer la surface la plus tendue qui passe par les quatre bords (qui interpole 
ces bords) : 

55 S(u,v)= (l-u)gl(v)+ug2{vW-v)n(u)+vf2(u)-[0-u)(A-v)P(0 } 0) 

+( 1 - u) vP(0 t 1 )+ u{ 1 - v) P(1 ,0)+ uvP{ 1,1)] ( 1 ) 
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Les sommets P(i, j) sont les quatre coins de la surface. 

On peut alors obtenir un maillage regulier NxM de cette surface par simple echantillonnage en prenant pour som- 
mets les points S(i, j) correspondant aux ui = i/(N-1), vj = j/(M-1), i=0 a N-1 , j=0 a M-1 . La formule 1 restreinteau points 
du maillage devient : 

sa,j) = d w (;>^ 92 (y >(i -jL )n {f)+J L m 
-[(i-^)(i-^)P(o,o )+ (i- 7 ^)^P(o,i) 

+ 7^< 1 -^i )P(1 ' 0)+ Afe/^ P(1 ' 1)] (2) 



Dans le cas illustre, les courbes fl, f2, g1, g2, sont obtenues a partir du logiciel GOCAD precite et sont done des 
lignes brisees. Pour obtenir un maillage regulier MxN, on echantillonne regulierement (en abscisse curviligne ou lon- 

15 gueur de segment) les courbes f 1 et f2 en N points, et les courbes g1 , g2 en M points, puis on applique la formule (2). 
On peut tout aussi bien obtenir un maillage a progression geometrique sur chaque direction, simplement en imposant 
cette contrainte au niveau du reechantillonnage des courbes. 

Ce maillage intermediate sauf cas exceptionnel ne convient pas car il n'y a pas de raison qu'il coincide a ce stade 
avec la surface qui a prealablement ete triangulee par le logiciel de surtacage. L'ajustement necessaire consiste a 

20 projeter les noeuds de la grille intermediate sur la surface triangulee orlhogonalement a celle-cL ou suivanl d'aulres 
directions , de fa9on a obtenir une nouvelle grille avec des noeuds places sur la surface triangulee. Quand les quatre 
courbes frontieres sont regulierement echantilionnees, la methode d'interpolation de Coons presente I'avantage de 
produire une grille reguliere en ce sens que sur chaque courbe isoparametrique tous les segments echantillonnes ont 
la meme longueur. Plus generalement la methode de Coons conserve les rapports de faille entre elements consecutifs 

25 dans une direction donnees. Ceci permet d'avoir par exemple des maillages reguliers ou respectant certaines pro- 
gressions dans la faille des mailles. Cependant la grille obtenue par cette methode connue perd sa regularity (ou le 
critere recherche) une fois qu'elle est projetee sur la surface triangulee. 

2) Rearrangement des noeuds de face 

30 

Pour retrouver le critere desire dans la distribution des noeuds de la grille projetee, une regularity satisfaisante 
par exemple, il faut deplacer les noeuds du maillage sur la surface triangulee. Pour ce faire, on procede a une relaxation- 
projection de la grille par ajustements successifs comme indique ci-apres (dans le cas ou le critere est la regularity), 
en relation avec les Fig. 1 , 2. 

35 Etant donne un bloc maille, le but cette relaxation est de rendre son maillage le plus regulier possible au sens ou 

n'importe quelle courbe isoparametrique du maillage soit regulierement echantillonnee. Pour ce faire, on deplace de 
facon iterative tout point du maillage de maniere a le positionner au mieux a egale distance de ses deux voisins sur 
chacune des trois directions (u,v,w) de voisinage. Soit P le point qu'on cherche a deplacer Pour chaque direction i, (i 
= u, v,w) (Fig!l7) : on trouve le point Pi qui se situe a egale distance des deux voisins de P sur la courbe de direction 

40 i, (i - u, v,w). Le point P est alors remplace par le barycentre des trots points Pi (i = u, v, w). 

Apres avoir deplace tous les points du maillage de la sorte, ceux-ci ne sont plus sur la surface triangulee. On les 
projette alors a nouveau sur la surface triangulee et on repete le processus de relaxation-projection jusqu'a atteindre 
la regularity souhaitee. 

Ceci est illustre sur I'exemple des Fig. 1, 2, correspondant a une face du cube deforme de la Fig. 3 qui est un 
45 exemple synthetique permettant de bien illustrer les differentes etapes de la methode. Cette surface a ete obtenue 
par simple deformation d'une face carree en preservant ses bords ( le bord de la surface reste un carre). Sur ces 
figures, on peut voir trois maillages 4X4. Le premier est un maillage plat issu de I'interpolation de Coons appliquee 
aux quatre bords (formanl un carre). Le deuxieme est le maillage fin donl les noeuds sont obtenus par projection de 
ceux du premier maillage. On constate que ce maillage a perdu de sa regularity. Le troisieme est le maillage regulier 
50 trace en gras et dont les noeuds sont situes sur la surface. II a ete obtenu a partir du deuxieme maillage apres une 
serie de projections et relaxations. La Fig. 4 montre le resultat de ces operations appliquees a I'ensemble des faces 
du cube deforme. 



3) Peuplement en noeuds de grille de chaque bloc dans son epaisseur : 

Pour prolonger en profondeur les grilles formees sur les faces des blocs (Fig.5), on utilise une generalisation au 
volume de la methode d'interpolation transfinie precedente. 

Etant donne 1 2 courbes parametriques f 1 (u), f 2(u ), f3(u), f4(u), gl (v), g2(v), g3(v), g4(v), hi (w), h2(w), h3(w), h4 
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(w) (0 2 u, v, w 2 1) definissant les 12 bords d'un bloc tet que montre dans la Fig. 16, une generalisation de la formule 
de Coons tournit une parametrisation volumique a Pinterieurdu bloc. La formule de Coons dans le cas volumique est 
la suivante : 

V(u,v,w) = ^-v)^-w)n{u)^0-v)wf2(u)+vwf3{u)+v0-w)f4(u) 
+( 1 - u)( '[ - w)g1 ( v)+( 1 - u) wg2( v)+ uwg3( v)+ u{ 1 - w) p4( v) 
+(1 - u)( 1 - v) M ( w)+( 1-u) vh2( w)+ uvh3( w)+ iv(1 - v) hA{ v) 

+ C(u,v,w) (3) 

avec 

C(l/,i/,w)=-2[(1 -u)(1 - v)(1 - iv)P(0,0,0)+(1 -v) wP(0,0,1 ) 
+(1-u)v(1-iv)P(0,1,0)+(1-u) i/tvp(0,1,1)+u(1-v)(1-iv)P(1 t 0,0) 
+ u(1 - v) wP(1 ,0, 1 )+ w( 1 - w) P(1 , 1 ,0)+ uvwPO ,1,1)) 

Les sommets P(i, j, k) ( i, j, k= 0 ou 1 sur la Fig 1 6) sont les 8 coins du bloc. 

En combinant les equations (1) et (3) on aboutit a une expression de la formule de Coons dans le cas volumique 
en fonction des 6 surfaces de bord du bloc, S1 (u, w), S2(u,w), S3(v, w), S4(v t w), S5(u,v), S6(u,v). Mais ici, les surfaces 
sont aussi supposees etre obtenues par interpolation de Coons a partir de leurs courbes de bord. La formule est la 
suivante : 

V( u, v, w) = ^[(1 - v)S1 ( u } w) + vS2(u, -u)S3{ v, w)+ uSA{ v, w) 

+(-\-w)SS(u,v)+wS6(u,v))+ 1 C(u,v,w) (4) 

Un echantitlonnage quelconque de V(u,v,w) fournit alors un maillage EF du bloc. 

Cette interpolation en profondeur donnerait une distribution reguliere de noeuds si elle etait appliquee a partir des 
noeuds des grilles formes sur les differentes faces du bloc. Ce n'est pas le cas et la distribution obtenue est genera- 
lement irreguliere car la methode d'interpolation en profondeur est appliquee aux noeuds reportes par relaxation sur 
les surfaces triangulees, comme decrit dans Tetape 2 precedente. 

4) Rearrangement des noeuds de bloc 

Aussi, pour retrouver une regularity suffisante (ou respecter un autre critere choisi), la methode selon ('invention 
comporte egalement une relaxation de la grille 3D de la maniere indiquee precedemment. Les Fig.5, 6 montrent le 
maillage interne du cube deforme avant et apres relaxation. Le maillage de la Fig. 6, on le note, possede une reguiarite 
bien meilleure que le maillage de la Fig.5. 

Quand chacun des blocs precedents de la structure etudiee a ete convenablement maille en volume, on reassem- 
ble tous ces macro-blocs apres leur reticulation 3D. 

La methode qui vient d'etre decrite peut etre utilisee pour maillerdes heterogeneites diverses telles qu'un chenal 
ou encore pour mailter radialement un puits devie. 

On a decrit des modes de realisation de ia methode appliquee au maillage d'une structure geologique. On ne 
sortirait pas du cadre de I'invention, cependant, en utilisant la methode pour le maillage de tout corps ou structure dont 
la geometrie est connue et imposee a priori. 

Revendications 

1. Methode pour generer un maillage 3D dans une structure dont la geometrie est connue et imposee a priori telle 
que par exemple une structure geologique definie a partir de donnees disponibles acquises par des mesures, par 
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des analyses et/ou des interpretations de donnees sismiques ou des diagraphies realisees dans des puits au 
travers de la structure, dans le but de realiser un modele representatif de cette structure, comportant une mode- 
lisation prealable des limites par des surfaces triangulees, enveloppant au mieux ces limites, et un f ractionnement 
de la zone en blocs ayant au plus six faces, dont les faces correspondent a des discontinues ou des failles 
principales de la zone prealablement triangulees, la methode etant caracterisee en ce qu'elle comporte : 

la subdivision de chaque bloc en un ensemble de cellules elementaires necessaires pour modeliser la zone 
par elements finis, par une reticulation ou maiilage en 3D de chacun des differents blocs obtenus, subdivision 
qui comprend : 

a) le maiilage de chacune des faces du bloc par interpolation a partir de ses bords, de facon a creer une 
grille intermediaire, sutvi d'une projection sur la surface triangulee des noeuds de cette grille intermediaire 
suivant une direction definie, de facon a generer des grilles de face; 

75 b) le peuplement de chaque bloc en noeuds internes par interpolation a partir des noeuds des diflerentes 

grilles formees sur les diflerentes faces triangulees de chaque bloc, de facon a obtenir une distribution 
reguliere des noeuds dans le volume de chaque bloc; et 

< c) un rearrangement dans le positionnement des noeuds des cellules elementaires creees, obtenu par 
20 relaxation dans le positionnement des noeuds, en fonction de crileres predefinis. 

2. Methode selon la revendication 1 , caracterisee en ce que Ton procede successivement a une premiere operation 
de relaxation dans le positionnement des noeuds sur chacune des differentes faces triangulees, le peuplement 
en noeuds internes etant effectue a partir des grilles relaxees, et ensuite une deuxieme operation de relaxation 

25 dans le positionnement de ces noeuds internes. 

3. Methode selon la revendication precedente, caracterisee en ce que Ton realise la premiere operation de relaxation 
par des phases iteratives de relaxation-projection. 

30 4. Methode selon I'une des revendications 1 ou 2 ou 3, caracterisee en ce Ton choisit comme direction de projection, 
une direction perpendiculaire a la grille intermediaire. 

5. Methode selon I'une des revendications precedentes, caracterisee en ce que Ton utilise un mode d'interpolation 
de type transfini. 

35 

6. Methode selon I'une des revendications precedentes, caracterisee en ce que Ton rearrange le positionnement des 
noeuds des cellules elementaires creees de facon a obtenir des mailiages reguliers. 

7. Methode selon I'une des revendications 1 a 5, caracterisee en ce que Ton rearrange le positionnement des noeuds 
40 des cellules elementaires creees de facon a obtenir un mode determine de progression de la taille des mailles. 

8. Methode selon I'une des revendications precedentes, caracterisee en ce que Ton procede a un reassemblage de 
tous les blocs. 
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